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Based on the notion of a construction process consisting of the stepwise addition of particles 
to the pure fluid, a discrete model for the apparent viscosity as well as for the maximum packing 
fraction of polydisperse suspensions of spherical, non-colloidal particles is derived. The model 
connects the approaches by Bruggeman and Farris and is valid for large size ratios of consecutive 
particle classes during the construction process. Furthermore, a new general form of the well-known 
Krieger-Dougherty equation allowing for the choice of a second-order Taylor coefficient for the 
volume fraction (^'^-coefficient) is proposed and then applied as a monodisperse reference equation 
in the course of polydisperse modeling. By applying the polydisperse viscosity model to two different 
particle size distributions (Rosin-Rammler and uniform distribution), the influence of polydispersity 
on the apparent viscosity is examined. The extension of the model to the case of small size ratios as 
well as to the inclusion of shear rate effects is left for future work. 



I. INTRODUCTION 



Multiphase flow is a very important but unexploited 
field of research according to the variety of unsolved ques- 
tions. In both nature and technology multiphase flow is 
rather the rule than the exception. The fleld of applica- 
tions includes spray, bubbly flows, process and environ- 
mental engineering, combustion, rheology of blood and 
suspension multiphase systems as well as electro- and mag- 
netorheological fluids. At this point of time, researchers 
neither agree on the cause of various effects nor on their 
theoretical description. In this paper, focus is put on sus- 
pensions. Depending on the length scale of observation as 
well as on the effects to be described, different approaches 
may be chosen. On the scale of individual particles in 
a microscopic approach, besides the strongly restricted 
possibilities for exact calculation [T] , SxOKESian dynamics 
([2, 3 ) and Lattice-BOLTZMANN methods ([HIS]) are em- 
ployed, for instance. These methods provide a means to 
investigate the mechanisms occurring in suspensions and 
to explain the origin of macroscopically observable effects. 
However, resolving individual particle requests a large 
computational effort and is therefore not applicable for 
engineering purposes. Increasing the observation length 
scale thus leads to reduced computational effort but also 
to loss of information because of coarser sampling. The 



EuLER-LAGRANGE-method uses groups of particles — so- 
called parcels — to represent the particle phase within the 
fluid carrier phase whereas the EuLER-EULER-method 
considers both phases as interacting continuous media 
(see e.g. [1]). Both methods are frequently used to in- 
vestigate transport, dispersion and reaction processes in 
dilute and dense suspension systems. In case a pure 
macroscopic description of the suspension is sufficient, 
one may model certain flow parameters of the suspension 
as a whole so that only, for instance, the volume frac- 
tion of the particle phase has to be determined during 
the computation to provide a basis for the calculation 
of macroscopic suspension properties such as the appar- 
ent viscosity. Here, micromechanics models ([71[S]) are 
well-suited, especially the notion of construction processes 
employed in the present work. 

In the following, we intend to describe the apparent vis- 
cosity as well as the maximum packing fraction of disperse 
systems by means of the volume fractions of the particle 
phase. Especially, we consider polydisperse suspensions as 
they are the most general case of dispersions. The appar- 
ent viscosity is a macroscopic quantity resulting from the 
presence of the particles. This quantity was first described 
by Einstein [T] for the dilute fimit, that is small volume 
fractions of the particle phase. Later, various attempts to 
extend the validity of the viscosity relation were made as 
reported in ([THH]). Taking into account that the apparent 
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viscosity is a macroscopic quantity, individual particles 
are not considered, rather the particle size is represented 
by particle size classes. The suspension is assumed to be 
built up from a finite number of such classes in the course 
of a construction process that allows for the calculation 
of the macroscopic suspension properties and thus needs 
to be described in detail. Since the suspension shows in 
general non-Newtonian behavior, the influence of the shear 
rate on the apparent viscosity will thus also be modeled 
in a future study. 

The present paper is organized as follows. In section 
some basics on apparent viscosity are provided along with 
an overview of recent work. Section IIIII is dedicated to 
a generalization of the viscosity correlation for monodis- 
perse suspensions. In section [TV] a construction process 
approach used to describe polydisperse suspension viscos- 
ity by monodisperse viscosity correlations is presented. 
Accordingly a model for the maximum packing fraction 
of polydisperse systems is outlined. The resulting model 
is then applied to two different particle size distributions 
demonstrating the effect of polydispersity on the apparent 
viscosity. Section |V] is devoted to conclusions. 



II. BASICS OF APPARENT VISCOSITY AND 
REVIEW OF RECENT WORKS 

Following the common quasi-single phase formulation 
we assume the rheological behavior of suspensions to be 
described by a stress tensor of the form (7] 

T = T(p,T,S) (1) 

where T and S denote the temperature and the symmetric 
part of the velocity gradient tensor S — [Vu -I- (Vu)-^] , 
respectively. We confine ourselves to the case of incom- 
pressible flow with a constant deformation history (so- 
called viscometric flow, see [TD]). Within the framework 
of continuum mechanics ( jlll 1121) these presumptions are 
met by the constitutive relation of the so-called general- 
ized NEWTONian fluid 

T = -pI + 2r,{T,j)S (2) 

with the generalized shear rate 

7 := \/2trS2 (3) 

which is equivalent to the second invariant of the symmet- 
ric part of the velocity gradient tensor S. The quantity 77 
expresses the apparent viscosity to be determined. 

Throughout the present work, we will disregard the 
existence of single particles but represent the particles 
summarily by the properties of so-called particle size 
classes. The presence of particles within the flow increases 
the viscous dissipation compared with the pure fluid phase 
with viscosity rjQ which leads to the observability of the 
apparent viscosity ?/. In order to isolate the influence of 



the particle phase on the apparent viscosity one defines 
the relative viscosity r]r as 

Vr - (4) 
Va 

The apparent viscosity is mainly dependent on the volume 
fraction 4> 

:= y ^7^'° (5) 

Vfluid + Vparticlo 

of the particle phase, where Vparticie and Vfluid denote 
the volumes of the particle and fluid phase, respectively. 
The volume fraction is sometimes called packing fraction 
as well. From experiments fl3J it is well-known that the 
relative viscosity monotonously increases with increasing 
volume fraction and exhibits a singular behavior at a 
value (j) < 1. The point where the relative viscosity 
diverges is commonly denoted as the maximum packing 
fraction (pc- By experiment, it has been found that for 
monodisperse suspensions (p^ coincides with the so-called 
random close packing state of the suspension which equals 

0c « 0.64 (6) 

for spheres (see |fj). A value for the maximum packing 
fraction of course also exists in the case of polydisperse 
systems but differs from the monodisperse case in general. 
We will propose a model for the polydisperse maximum 
packing fraction in section |IVE[ 

Suspensions having a fixed particle size distribution 
but a variable volume fraction of the particle phase be- 
have quantitatively similar to the monodisperse suspen- 
sion described above. The quantitative changes concern 
the maximum packing fraction 0c and the slope of the 
curve rjr — rjri'P)- These effects have been examined sys- 
tematically in |14j . see Figure [l] Therein the relative 




volume fraction <j> 



Figure 1. Dependence of the relative viscosity on the volume 
fraction and the width of the particle size distribution from 
the experiment |14) : narrow (ND), moderately broad (BDl) 
and very broad (BD2) distributions 

viscosity is depicted as a function of the volume fraction 
for three distributions of various width. Two trends can be 
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observed. On the one hand, the viscosity of a polydisperse 
system with a broad size distribution always hes below 
the viscosity of a system with a narrow distribution, thus 
below the viscosity of a monodisperse system as well (see 
also [T3] [121 US]). On the other hand, the maximum pack- 
ing fraction increases with increasing polydispersity, that 
is with the broadening of the particle size distribution. 

As already mentioned, the fundamental work on the 
apparent viscosity of disperse systems has been written by 
Einstein (see [T], correction [H]). Therein, the Stokes 
equation 



r]Au = — Vp 
P 



(7) 



is solved in a three-dimensional dilatational flow around 
a spherical particle at rest. Afterwards, the solution 
is transferred to the case of a suspension with a finite 
number of particles (volume fraction (p). The dissipation 
change due to the presence of the particles leads to the 
well-known Einstein relation 



TJr ^1 + 2.5(1) 



(8) 



Though fluid inertia is neglected in the Stokes equa- 
tion (|7]), the equation ^ serves a an exact limit for dilute 
suspensions, that is for (f> 0. The Einstein relation is 
commonly considered as the first-order TAYLOR series ex- 
pansion of every correlation between viscosity and volume 
fraction [7]. We denote the (/>— coefficient as first-order 
intrinsic viscosity [771] and analogously [77^] as mth-order 
intrinsic viscosity, so that we can write down the TAYLOR 
series expansion of the relative viscosity in the form 



T]r-1+ [m] <P + [V2] (f>^ + --- + [Vn] 0" + 

oo 

= 1 + ^ [ry„] 0™ 

m— 1 



(9) 



This representation will be used later in this work. In the 
literature a great number of equations of the form 



(10) 



are provided. Tables [T] and [IT] summarize the most impor- 
tant correlations considering neither their derivation nor 
the underlying models. Most of the equations are empiri- 
cal or semi-empirical, only few are exact. The reason for 
the existence of such a great number of different correla- 



tions is that relations of the form ( 10 1 do not cover the 



entire parameter space governing the physical problem. 
In [7J this is expressed by the formulation 

r]r — r]r{(j), all other details of microstructure) (11) 

Clearly the parameter space must be confined to allow 
for useful modeling and thus the range of validity has to 
be confined a priori. In the present context, the correla- 
tions (11), denoted as viscosity relations in the following. 



a. Series expansions with respect to the volume frac- 
tion for cf) <^ 1 according to the equation 

While the Einstein relation ^ with the intrinsic vis- 
cosity [rji] = 2.5 is commonly accepted as the first order 
series expansion of the relative viscosity r]r of suspensions 
with spherical particles, there is no unique value of [772] 
because of a strong case-sensitivity of this parameter. In 
Table [I] some examples taken from the literature are listed. 

b. Correlations for the entire range of (p: 

The correlations listed in Table |ll] are intended to be 
valid over the entire range of 0. They all coincide with 
respect to a singular behavior at the point (p — (pc, that is 
when the maximum packing fraction is reached. For ^ — ^ 
the equations by Frankel-Acrivos and QUEMADA both 
do not reduce to the Einstein relation. Therefore, in [7] a 
correction to the Frankel-Acrivos relation is proposed 
in the form 



?7r = 1 + 



(i) 



1/3 



20 



1/3 



(18) 



9 



1 



2/3- 



Figure [2] shows the relations Q to ([T7| from Table |lT] 
for the parameters cpc = 0.64 and [771] — 2.5. In the next 
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Figure 2. Correlations between relative viscosity rjr and volume 
fraction p for < p < pc according to the equations (12 I 
to (Il7| with pc = 0.64 and [r?i] = 2.5 



can be classified in two groups: 



section, an attempt to generalize the viscosity correlation 
for monodisperse suspensions is presented. 
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Author(s) 



Reference Annotations 



Batchelor & Green (1972) [S] 

Batchelor (1977) [9] 

Bedeaux et al. (1977) [9] 

CiCHOCKi & Felderhof (1991) [5] 

CiCHOCKi & Felderhof (1991) [9] 

Kim & Karrila (1991) [H] 

Krieger & Dougherty (1959) [19] 

MooNEY (1951) [7] 



BROWNian motion neglected, random spatial particle 
distribution 

BROWNian motion included, random spatial particle 
distribution 

formalism in wave number space 

BROWNian motion neglected, Smoluchowski equations 
BROWNian motion included, Smoluchowski equations 
STOKESian dynamics similar to Batchelor (1972) 
second-order Taylor coefficient of the 
Krieger-Dougherty relation ([T2| for (pc = 0.64 



second-order Taylor coefficient of the Mooney relation ( 13 \ for (j>c = 0.64 7.03 



5.2 

6.17 

4.8 

5.00 

5.91 

6.95 

5.08 



Table 1. Values of the second-order intrinsic viscosity [7)2] taken from the literature 



Author(s) 



Reference Equation 



No. 



Krieger & Dougherty (1959) [19] 



= 1- 



-lvi]4>c 



Mooney (1951) [16] rir 

Frankel & AcRivos (1967) [H] rjr 

Filers (1941) [7] rjr 

Quemada (1977) [ID] ijr 

Robinson (1949) [7] ijr 



(12) 
(13) 

(14) 



-P(l*) 

[i + |M(i47^;)]' (15) 

i^-J^y (16) 

l + ['7i](i4z^) (17) 



Table II. Correlations between relative viscosity 77,. and volume fraction <j) for < < 



III. GENERALIZATION OF THE VISCOSITY 
CORRELATION FOR MONODISPERSE 
SUSPENSIONS 

As shown in the previous section, there are two main 
types of viscosity correlations, namely polynomial and 
closed correlations. Polynomial correlations are well 
suited for describing the low-concentration range but 
do not show divergence for <j> ^ (f>c. Closed correlations 
diverge for (p ^ (pc and therefore seem to be suited for 
the whole concentration range, but cannot show proper 
asymptotic behavior for </() — > because the second order 
TAYLOR-coeffrcient [772] is determined a priori through the 
viscosity relation. So in the following section we combine 
the low-concentration behavior of polynomial correlations 
with the high-concentration behavior of a special closed 
viscosity correlation based on the well-known Krieger- 
DOUGHERTY equation. 



A. Derivation of the viscosity correlation 

We derive a correlation that in addition to the prop- 
erties of the closed correlations in Table |ll] allows for 
choosing the second-order Taylor coefficient [r/2] . This 



value depends on particle interactions and may differ con- 
siderably from case to case. So the modified viscosity 
correlation shall include the three parameters [r]i], [r/2] 
and (hr- 



Starting from the Krieger-Dougherty equation ( 12 1 



-[»?l]0c 



(19) 



we assume the modified correlation to be of the form 

-C 



■qr = A + B{l^j^ 



(20) 



This ansatz has to fulfill the conditions 



r]r{(j) ~ 0) = 1 and 
?7r(</)^0) = l-F[ryi], 



Condition (21) implies that 



A + B 



(21) 
(22) 



(23) 
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The second condition ( 22 1 is equivalent to 

= [m] 



dr]r 



6=0 



= 2 [,72] 



(24) 



(25) 



(t>=0 



In order to ensure the fulfillment of this condition we 



calculate the first two derivatives of the ansatz ( 20 1 with 



respect to cj) and additionally use the equation |23|. We 
obtain 



df]r 



BC 



dcf) (j>c 



1 - — 



-C-l 



(C + 1) 



BC 



-C-2 



(26) 



(27) 



and evaluate the derivatives at the point ^ = 0. It follows 
from the conditions ( 24 1 and ( 25 ) 



A=l 



B = 



: [m 



2<f>c [m] - [m] 



^^c [V2] - [Vl] 

2cj)c [V2] - [Vi] 



(28) 
(29) 



.[Vi] 



so the final form of the modified viscosity equation ac- 



cording to the ansatz ( 20 1 is 



?7r(0) = 1 - 



2(/'c [m] - bii 



24^c[V2]-{vi] 



(31) 



B. Special cases of the viscosity correlation 



Since equation (31 1 has been derived under rather for- 



mal considerations than from a physical point of view, it 
has to be understood as an empirical equation to be fitted 



to experimental data. Because equation (31) reduces to 



the Krieger-Dougherty equation al9« for the choice 



(1 + [m] cf>c) [vi] 

20c 



(32) 



it is of course superior to the latter for fitting purposes. 
As a consequence, we do not learn anything from fitting 
equation (31) to experimental results for testing. 

Instead we examine if the closed correlations in Table HIl 
can be reproduced by the modified viscosity equation. 
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[1)2] = 7.9, cfc = 0.54 



MooNEY relation (13) 
(Ac = 0.64 



Equation (31) 

[1)2] = 4.4, 0e = 0.64 



Robinson relation (17) 
0c = 0.64 



-1_ 



0.1 



0.2 0.3 0.4 

volume fraction (t 




Figure 3. Comparison of equation (31 1 with the relations by 
Mooney (13 1 and Robinson (171; [771J = 2.5 



From Figure |2] we chose two extreme examples, namely 
the relations by MoONEY and ROBINSON (equations (13) 
and (17)). Figure [s] shows the results. We find gooc 
agreement in both cases. However, it was necessary to 
decrease the value of the maximum packing fraction (j)c 
in order to reproduce the behavior of the Mooney equa- 
tion because of its exponential variation compared with 
(30) the algebraic variation of the equation (31). The remain- 



ing correlations in Table [11] vary algebraically and may 
therefore be reproduced by the equation (31 ) with high 
accuracy. 

Because of these reproduction properties we will in the 
following consider the modified viscosity equation (31 1 as a 



general viscosity correlation for monodisperse suspensions. 



IV. DEVELOPMENT OF A VISCOSITY 
CORRELATION FOR POLYDISPERSE 
SUSPENSIONS 



In this section we develop the polydisperse viscosity 
model based on the notion of a construction process. 
This approach is first exactly described in section |IVC| 
Subsequently in section |IVD| the construction process is 
applied to the determination of relative viscosity. Then, a 
model for the maximum packing fraction of polydisperse 
suspensions is developed in section |IV E| to complete the 
viscosity model. A graphical scheme provided in Figure [6] 
may serve for the reader's guidance during the calculation. 

It is important to note that the basic model developed 
in this section can only be applied to polydisperse suspen- 
sions with large diameter ratios of consecutive particle size 
classes in the construction process. Work on an extension 
of the model to small particle size ratios is in progress. 
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Starting point: The differential Bruggeman 
model 



The diflFerential Bruggeman model (see also |2T] and 
more detailed |[8j) makes it possible to derive a closed 
viscosity relation for the full concentration range starting 
from the EiNSTEiN relation. The Bruggeman model 
is also known as Differential Effective Medium approach 
(DEM). A generalization of the DEM approach is pre- 
sented in [22 . The Bruggeman model is based on the 
notion that an infinitesimal volume fraction of particles (p* 
is added to an existing suspension with effective viscos- 
ity 77 and volume fraction </>. In the course of this addition 
it is assumed that the existing suspension can be treated 
as a homogeneous medium. This can only be valid if the 
newly added particles have a large diameter compared 
with the particles already present in the suspension. 

We now ask for the change in effective viscosity due to 
the newly added volume fraction (j)* . The infinitesimal 
volume fraction of the newly added particles in the re- 
sulting suspension is, according to |21| (compare the later 
relations ((45| and ([48|) 



1 + . 



1 - 



(33) 



Because of the small size of the volume fraction (pp we 
may use the Einstein relation to describe the change in 
effective viscosity by 



77 -I- d?7 ?7(1-|- [771] 4)p) 



(34) 



By inserting the equation ( 33 1 into the relation ( 34 ) we 
obtain 



rj + drj = rj I 1 + [rj 



d0 
1^ 

Simplification and separation of variables yield 

dr/ d(j) 
[771] 77 1-0 



(35) 



(36) 



We integrate the equation ( 36 1 under the initial condition 
77(0 = 0) = 7^0 for the pure fluid and find 



V '70 (1 - (/)) 



-[vi 



(37) 



The equation (37 1 is known as ROSCOE equation. Since 



the Bruggeman model requires a large diameter ratio of 
consecutively added particle classes, the suspension must 
therefore consist of a solid phase that can be divided 
into particle size classes of large diameter ratios. This 
structure is called hierarchical, see also [22j and [81 So 
the differential Bruggeman model can be applied to 
poly disperse Suspensions. 

Another important assumption of the differential 
Bruggeman model is the vahdity of equation (l34|. The 



volume fraction of newly added spheres in equation ( 33 1 
has to be small enough for the Einstein relation to be 
valid. 



We note that the volume fractions (j)* and accord- 
ingly (f)p may be finite in principle. However, by intro- 
ducing the infinitesimal increment d77 into the differential 
Bruggeman model the volume fractions (p* and pp are 
required to be infinitesimal. The advantage of this limita- 
tion is the possibility to derive the closed equation ( 37 ) . 



B. Assumptions 

In the following we assume the suspension to consist of 
a solid particle phase suspended in a fluid carrier phase 
with constant viscosity 770. The particle phase consists of 
spheres with different diameters di that can be categorized 
in a finite number 77, of size or diameter classes. So with i 
as the index of the size class we have i — I . . .n. The size 
classes shall be sorted by diameter in ascending order, so 
that di < di+i. The ratio of two consecutive diameters 



di+i 
di 



(38) 



should be larger than 7 according to so that the ex- 
isting suspension behaves like a homogeneous medium 
towards the newly added spheres. In the completed sus- 
pension resulting from the construction process the ith 
size class occupies a volume Vi while the fluid phase oc- 
cupies the volume Vf. So the total volume Vp of the 
suspension is given by 



Vt = Vf 



E 



Vrr, 



(39) 



We assume that the suspension can be entirely described 
by volume fractions, so that no higher moments of the 
so-called indicator function (see [8J) have to be considered. 
By the way, the indicator function is defined on the entire 
space occupied by the suspension and is unity at all points 
belonging to phase 1 (e.g. the solid phase) and zero for 
all of the remaining points (e.g. the fluid phase). The 
function is commonly used for a statistical description 
of phase interactions by means of its moments. The 
lowest-order moment is given by the particle phase volume 
fraction. Confining ourselves to a description based on 
volume fractions corresponds to the assumption of an 
isotropic and homogeneous suspension throughout the 
control volume. 

As already outlined in section ^ the total volume 
fraction is defined by 



V 



particle 



Va 



uid 



V 



particle 



V" V 
Vt 



(40) 



Analogously it is useful to define volume fractions of single 
size classes, both during the construction process and in 
the completed suspension. 
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C. Construction process 

The models for the apparent viscosity and the maxi- 
mum packing fraction that are developed in the following 
sections are based on the notion that the suspension is 
constructed by successive addition of new size classes. We 
call this process the construction process. In the following 
we generalize the considerations in [7] and especially [25]. 
There are two possible approaches for the construction 
process shown in Figure |4] 

Variable total volume: In this case the volume of the 
fluid phase V/ is held constant during the construc- 
tion process, so that the total volume of the suspen- 
sion increases with each step until the suspension 
occupies the final volume Vt- The construction pro- 
cess thus only consists of additions of size classes. 

Constant total volume: In order to keep the total vol- 
ume of the suspension constant throughout the con- 
struction process, it is necessary to extract a suspen- 
sion volume with a size equal to the added particle 
volume in each construction step (see also [H]). So 
the extracted volume represents the composition of 
the existing suspension. 

Though we will show in a later section that both ap- 
proaches are equivalent, we will choose the constant- 
volume approach for reasons of simplicity. The homog- 



Figure 4. Scheme of the construction process for a bidisperse 
Suspension (n = 2) at variable (above) and constant volume 
(below) including homogenization 

enization step in Figure |4] symbolizes the fact that the 
existing suspension acts as a homogeneous medium to- 
wards the newly added particles because of the large 
diameter ratio. 



1. Construction process at variable total volume 

In the case of variable total suspension volume the 
added volumes of the respective size classes are identical to 
the volumes of the size classes in the completed suspension 
at the end of the construction process. In contrast, this 
is not the case when the total volume is constant. 

First we show explicitly the addition of the first two 
size classes with volumes Vi and V2 to the fiuid phase 
with volume Vf . The total volume fraction of the particle 



phase after the ith construction step (pi is given by (see 

e.g. m) 







= 



Vf+0 

Vi 
Vf + Vi 

V1 + V2 

Vf + Vi + V2 



(41) 



A general [i -f l)th construction step can thus be written 
in the form 



V-i-Hl y 



— >• 



i+1 



(42) 



If we relate the added volume Vi+i to the total suspension 
volume before the addition of Vi+i and call this ratio (pl+i 



(43) 



we can modify the expression for the total volume fraction 
of the particle phase (f)i+i in equation (42) so that 



Ui+l 



1 



(44) 



The concentration of the (z + l)th size class after the 
{i + l)th construction step is consequently given by 



V, 



(45) 



The volume fraction (?!)p.i+i in equation (45 1 will be impor 



tant for the calculation of the apparent viscosity because 
it serves as an argument in the viscosity relation (compare 
equation (34l). 

For a different formulation we denote the change in 
total volume fraction during the (i + l)th construction 
step as 



Using equation ( 44 1 we first obtain 



1 



and finally under consideration of definition ( 45 ) 



1 



(46) 



(47) 



(48) 



2. Construction process at constant total volume 

If the total suspension volume is held constant through- 
out the construction process, in every construction step 
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we have to extract as much volume as is added with the 
new size class. The extracted volume represents, as al- 
ready explained on page [7] the composition of the existing 
suspension and therefore contains particles of smaller size 
classes (see also [22J). This implies that the volume of the 
respective size class added in each construction step has 
to be larger than the volume of this size class contained 
in the completed suspension. We notice that the last 
size class (the largest particles) is an exception and is 
therefore added with its final volume. 

The volume added in the ith construction step is thus 
called V* in order to be distinguished from the volume Vi 
of the ith size class in the completed suspension. The 
constant total volume is still called Vt- 

Now we consider the construction process in detail. The 
first addition implies a simple change in the total volume 
fraction of the particle phase (jjf. 







= 



V{ 



(49) 



Vt Vt 
The more complex second construction step is given by 



Vt 



V{ - V, 



2 Vt 



+ v,* 



Vt 



(50) 



In the equation ( 50 ) the term in brackets represents the 



volume of the first size class still present after the second 
construction step. Therein the term V2V1/VT describes 
the loss of volume of the first size class due to the necessary 
extraction of volume. If we introduce the notation 



Vt 



(51) 



we may analogously to the equation ( 45 1 express the 
concentration of the {i + l)th size class in the existing 
suspension by 



(52) 



Unlike in the case of variable total volume we obviously 
do not need to distinguish between the volume frac- 



tions 



and 



in the equation ( 50 ) we find 



By rearranging of the expressions 



and 




(53) 
(54) 



In the equations ( 53 1 and ( 54 1 we have introduced a nota- 



tion for the volume fractions of the individual size classes 
during the construction process. The representation A0^. 
refers to the volume fraction of the fcth size class after 
the ith construction step, that is after the addition of 
the ith size class. So i means an index and no expo- 
nent. This should cause no confusion because the volume 



fraction will always occur linearly in all of the following 
expressions. 

It can easily be shown, using the equations (53 
and ( 54 1 , that the total volume fraction after the thirc 



construction step is given by 

<t)3 = 02(1 - ^3) + <f>*3 

= 0J(1 - cf,;){l - 0S) + </.;(! - 05) + 03 



(55) 



The underlying pattern can be recognized clearly and so 
we may generalize intuitively: 



t+l = 0l(l - 01+1) + ' 



(56) 



Table [111] schematically outlines the construction process. 
By introducing the difference 50i+i from the equation (46 1 
into the equation ( 56 I , we find 



h + S(j)i+i = 0^(1 - 0*+i) + 0-+1 sowie 



')P.l+l = 0*+i = 



1 



(57) 
(58) 



which is formally identical with the equation ( 48 1 . At 



tention should be paid to the fact that the volume frac- 
tions 0*^2 used in the cases of variable and constant 
volume, respectively, are not identical. The proof of the 



equivalence of the equations ( 48 1 and ( 58 1 is provided 
in appendix |A] for clarity. From this proof we draw the 
conclusion that the construction processes at variable and 
constant total volume are equivalent. The only difference 
within the framework of our dimensionless representa- 
tion, which is a representation by volume fractions, lies 
within the equation (A6). We choose the case of con- 



stant total volume because of the intuitive meaning of 
volume fractions originating from the constant volume Vt- 
Considerations involving volumes may therefore easily 
be transferred to the notion of volume fractions. In the 
case of variable total volume the total volume contrarily 
depends on the amount of added volumes. 

Representation of the quantities by the volume 
fractions A0fc 

Up to now, the equations contain quantities influenced 
by the construction process. However, in the case of 
calculation we only know the volume fractions of the 
size classes in the completed suspension. It is therefore 
useful to express the quantities describing the construction 
process by the composition of the completed suspension. 
For simplification we define 



A0fc A0" 



(59) 



These volume fractions in the completed suspension are 
given by 



A0fe = 



Vk 



Vf+El^.Vm 



Vt 



(60) 
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After Step 


Volume Fraction of Particle 
Phase 


SC 1 
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02 


Aoi? 






3 


03 




A<j!.i A<j!.i 




n 


4>n 


A0r 


A0? A0? . 


. AC 



Table III. Scheme for the volume fractions during the construction process (SC = size class) 



From the equations ( 53 1 to ( 55 1 we deduce a general 
expression for A0^: 



^'^k=<Pl n (l-'^™) ffirfc<z (61) 

rn—k-\-l 

A(/)fc fiir fc = i (62) 



We recall that there are no volume fractions with k > 



(compare Table III). If we succeed in calculating the 
volume fractions we are in a position to calculate 
the volume fractions of the individual size classes A(^^ 
and the total volume fractions after each step 4>i from the 



equations (61 1 and (62 1. Per definition we have 



k=l 



(63) 



(consider the equation ( 55 1 as an example) . Writing 



down the volume fractions of the last size classes in the 
completed suspension and using the definition (p9| as well 



as the equations (61) and (62) reveals the possibility to 
calculate the volume fractions 0^: 

A0„ = (f>n 



A^n-i 
A0„_2 = 
A^n-S = . . . 

So we find by recursive insertion that 

(f>n = ^(l>n 
,* A(/)„_i 



(64) 
(65) 
(66) 



(67) 



1 



A</>„ 

A0„_2 



1 ~ A(?!)„ - A0„^i 



(69) 



This may be generalized in the form 



We now introduce the abbreviation 

71 



(70) 



(71) 



into the equation ( 70 I and find 

1 jl* Z-^m—k 



A0„ 



1 - E™=fe+i A'?^" 

Therefore, we may write 

A0*, = </)*(l-0*+i)(l-0*+2)...(l -</.*) 

,* ffe+l crfc+2 CTfc+a cr.j 

= <Pfe 

Crk+2 CTfe+S Crfe+4 CTi+l 

,* CTfe + l 

= (/-fc- — 



(72) 



(73) 



instead of the expression (61). Using the result (70) and 



the definition (71) it follows from the equation (73) that 



A0l = 



A0fc 

1 - El=,;+i ^(l^n 



(74) 



Combination of the equations (|63| und (|74|) finally yields 

(75) 



1 - Em=i+1 

So we have represented all of the quantities occurring in 
the construction process by the volume fractions in the 
complete suspension. 



D. A discrete model for the relative viscosity 

As we have already noted, the differential Bruggeman 
model lacks any information about the volume fraction 
of the individual particle size classes. For that reason we 
have described the construction process of the suspension 
in a discrete form in section |IVC[ In preparation for the 
development of the discrete viscosity model we need to 
make a connection between the Bruggeman model and 
the maximum packing fraction. 



1. Introduction of the maximum packing fraction into the 
differential Bruggeman model 



The ROSCOE equation (37) diverges as the total vol 



ume fraction approaches unity. In a real suspension the 
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achievable value of (j) is limited by the maximum packung 
fraction. In order to introduce the maximum packing 
fraction into the differential Bruggeman model we pro- 
ceed in a way proposed in [21 . A quite similar way can 
be found in |i9j . In both publications it emerges that the 
notion of maximum packing fraction is introduced under 
little convincing considerations. 

The approach followed in [21; consists of modifying 
equation ( 33 1 by using the maximum packing fraction (f>c 
in the description of the volume fraction 0p. Therefore, 
it is supposed without derivation that 0p may be written 
as 



(76) 



In combination with the equation ( 34 1 one finds 



ri + di] = rj \ I + [r)i 



d0 



1 



(77) 



Rearrangement of the equation ( 77 1 analogously to the 
equation (36) yields 



dr] 



(78) 



which can be integrated under the initial condition ry((/) — 
0) — 7/0 for the pure fluid. This results in the Krieger- 
DOUGHERTY relation ([12 1 



-[')l]'/'c 



(79) 



In [5] the introduction of a so-called effective volume 
fraction leads to the equation 



-[vi 



(80) 



The expressions ( 79 1 and ( 80 1 obviously differ with respect 



to their exponents only. Both results show divergence 
for — 1 but lack a physical rationale for the respective 
approaches. Therefore, we will state a different principle 
in the next section. 



2. Introduction of the maximum packing fraction into the 
discrete construction process 

It would formally be possible to transfer the modifi- 
cation (76) to the discrete construction process, that is 



the volume fractions ( 70 1 , and employ the result for the 
viscosity calculation. In the following it will be explained 
why this approach cannot be valid in general. Partially 
anticipating the later viscosity calculation, we raise the 
two following points. 



Firstly, the construction process described in sec- 
tion IV C| is by no means dependent on the particle geom- 
etry. This is emphasized by the notion of homogenization 
between two construction steps. In contrast, the maxi- 
mum packing fraction is strongly influenced by the particle 
geometry. So it would be artificial to introduce this quan- 
tity into the description of volume fractions during the 
construction process. 

Secondly, the consideration of the volume fractions 
during the construction process is independent of the 
physical quantity that is calculated (here: the viscosity). 
It does not make any difference whether one calculates the 
viscosity or, for instance, the electric conductivity (or both 
at the same time). In both cases the construction process 
is constituted by the same volume fractions. Only through 
the employed relation between the volume fractions and 
the change in the physical quantity of interest parameters 
like the maximum packing fraction are included. This 
will be the approach followed during the later viscosity 
calculation. 

The above considerations imply that the differential 
Bruggeman model only allows for the derivation of the 
ROSCOE equation ( 37 1 because as a consequence of this dif- 



ferential approach the right-hand side of the equation ( 35 1 
may only consist of a linear expression (the Einstein rela- 
tion) that cannot contain the maximum packing fraction. 
So the approaches presented in [51] and [5] are formally 
possible but physically questionable. 



Notation for the maximum packing fraction in 
the construction process 

At this point it is necessary to introduce a distinct 
notation for the maximum packing fraction in order to 
avoid misinterpretations. The calculation of the maximum 
packing fraction will be conducted in section p^V E[ We 
have to distinguish between (jjc as a parameter within the 
correlations listed in Table |IT] and the modeled maximum 
packing fraction that may change during the construction 
process as well as under the influence of the shear rate. 
We choose the following notations referring to [23j: 
'^Tk '■ denote as ipl^f, the maximum packing fraction 
of a polydisperse suspension consisting of k size 
classes after the ith construction step (ith line in 
Table [ml . 

We write ipc '■= fxi fo'^ monodisperse packing 
fraction. The monodisperse maximum packing 



fraction ip^ is constant throughout the construc- 
tion process and thus carries no upper index i. 



This notation is visualized in Table IV (compare Table III ) 



After each construction step the maximum packing frac- 
tion is newly calculated according to the new composition 
of the suspension. This is conducted within a recursive 
process (arrows in Table IV I. Starting from the monodis- 
perse maximum packing fraction tpc all previously added 
size classes are taken into account and so for each step 
the value t/?^j^i is calculated. We will need this value for 
the equation (88) still to derive. 
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Table IV. Notation for the maximum packing fraction during the construction process (SC = size class); the arrows indicate 
that the far right values are calculated recursively from the previous values (for details of the calculation see section IV E I 



3. Transition to a discrete relative viscosity model 

We now build on the description of the discrete con- 
struction process given in section |IVC| and focus on the 
effect of the construction process on the relative viscosity. 
When during the construction process the {i + l)th size 
class is added, the total volume fraction of the particle 
phase changes from (pi to (f>i+i. This is associated with a 
change in apparent viscosity from rji to rji+i- To empha- 
size the analogy to the differential Bruggeman model, 
we temporarily confine ourselves to the linear Einstein 
relation as a description of the viscosity change and sub- 
sequently in section |IVD 3| we will extend the model to 
higher orders. 

According to sections |IVC 2| and |IV A| at constant total 
volume the volume fraction of the newly added (i + l)th 
size class related to the resulting suspension volume is 



So the new apparent viscosity is given by 



Vi (1 + [m] <Pp,t+i) 



(81) 



(82) 



In the equation ( |82| we have used the volume frac- 
tion 0p,i+i related to the suspension volume after the 
{i + l)th construction step according to the notion of the 
construction process. The volume fraction 4>*^i may not 
be used in general because the volume fractions (j)p^i+i 
and (f>*_^i only coincide in the case of constant total sus- 
pension volume (see section IV C 2 1 . All the following 
relations containing 0p,i+i may thus be used both at 
variable and constant total volume. 



In the latter case, the equation ( 82 1 can be cast in a 



different form employing the volume fractions A^j, from 
the equation ( 70 1 : 



77,+i = 77, 1 + [771] 



i+i 



(83) 



This equation is a recursive formula that can alternatively 
be written in an explicit form. However, we note that a 
numerical evaluation has to be conducted in a recursive 
way because the explicit formula offers no simplification 



compared with the recursive formula. The explicit form 
of the equation ( 83 I is given by 



i+l 



mY[ 1 + [m] 



A0„ 



(84) 



with (pQ = 0. For i + 1 — n we find the apparent viscos- 
ity r] — Tjn of the complete suspension (it is well known 
that X]ILri-(-i 2; = for all x, which may be important for 
a numerical implementation). 



Extension of the discrete model to higher order 
terms oi cpp 



The differential equation (78) contains only the first 



order Einstein relation because of the infinitesimal char- 
acter of dcj). In the case of the discrete model represented 



by equation ( 84 ) it is not necessary to confine oneself to 



linear terms. Therefore, we are allowed to describe the 
modification of the apparent viscosity more accurately 
by higher order terms of (pp^k- So the equation (84 1 can — 
according to the general expression 



Vi+i = Vq 



771—1 

— be extended to 



(85) 




(86) 



which according to the equation ( 83 1 may alternatively 
be written as 



Vi+i = Vi 



1 + iVi] 



A(j)„ 



■ + 



(87) 
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Comparison with relations present in the liter- 
ature 

In the following we draw a connection between the equa- 
tion (86) and relations existing in the literature. Selecting 



the coefficients from the Taylor expansion of a viscosity 
relation ((31), for instance) at the point ip ~ for the 
coefficients [ry^.] in the equation (85) and formally consid- 



ering an infinite number of terms, we may substitute the 



bracketed term in the equation ( 85 1 with the viscosity 
relation itself. This step requires the introduction of the 



maximum packing fraction. Using the equation (31) we 
find 



i+i 



m— 1 



1-1 



2^t™ [V2] - [Vi 



which only consider a relative but not an absolute influ- 
ence of the particle radius. Therefore, we will not use this 
approach in the current work. 



In later sections we will use the result ( 90 1 as a general 



formula for the apparent viscosity of hierarchical polydis- 
perse suspensions (compare section IV C 1. It is important 
to note that the extension of the equation ( 86 1 to ([90| 



has been conducted by means of considering the complete 
Taylor expansion of relation (31), suggesting the ad- 



missibility of arbitrary volume fractions of the individual 
size classes which is uncertain regarding the notion of the 
construction process. Since it is not possible to state a 
limiting value of 0p,i+i for the validity of the result (90), 



we will though use this equation at least as a reasonable 
approximation also for higher values of 4>ps+i- 

To sum up, we state that the equation (86) resulting 



form the construction process is at least justified for the 
first two orders of (pp^i+i, while the Farris model (cor- 
responding to the equation ( |90[ )) provides a useful but 
possibly unjustified extrapolation for high values of (j)p^i^i. 



E. Determination of the maLximum packing fraction 



In the equation (|88| refers to the maximum pack- 

ing fraction after the mth construction step that will 
be calculated in section [IV E This notation has already 



been outlined in section IVD2[ T he general recursive 
representation of the equation ( 88 ) is given by 



PA+l 



(89) 



and more specifically using the viscosity relation (31 ) 




(90) 



Despite the use of the variable maximum packing fraction 
the equation ( 88 1 corresponds to the so-called Farris 
model [24| also referred to in [23J and [25j. So we have 
found an interesting connection between the differential 
Bruggeman model and the Farris model which is drawn 
by the discrete construction process. Remarkably, the 
approach followed in |2S] includes a viscosity relation 
which depends on the particle size class, that is the par- 
ticle radius. Although this additional degree of freedom 
is useful for fitting purposes, it is inconsistent with the 
assumptions of the models by Bruggeman and Farris, 



In the following we will present two models for comput- 
ing the polydisperse maximum packing fraction. The first 
model has been proposed by FURNAS (see also [53] )■ Al- 
though this model is insufficient for reasons to be given in 
section |IVE 1 we will use the underlying considerations 



for developing an improved model in section |IVE2[ The 
approaches differ with respect to the underlying treatment 
of the particle size distribution. 

During the calculation we will use the quantities occur- 
ring in the construction process according to the tables III 
and |IV[ For instance, after the third construction step the 
values Lp^2 ^iid tp^^ are being calculated starting from 
the monodisperse maximum packing fraction ip^^^ — ipc 
and considering the existing volume fractions A0f , A(/)2 
and A0|. 



According to section 



IVD2 



we denote as 



the maxi- 



mum packing fraction after the ith construction step so far 
containing k size classes while the completed suspension 
consists of n size classes. 



1. The Furnas model 

The Furnas model is based on the assumption that 
the state of maximum packing fraction is constructed 
successively from size classes with decreasing diameter. 
At first, the larger spheres fill the entire available suspen- 
sion volume with the monodisperse maximum packing 
fraction ipc- So the total volume fraction is 



^1 = Vc 



(91) 



Subsequently, the remaining volume fraction (1 — 
is filled by the spheres with the next smaller diameter. 
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also to a realizable part of ip^- So the additional volume 
fraction of small spheres ^p\x given by 



we conclude that 



= (^c + (1 - (/Jc)Vc = 1 - (1 - 'Vcf 



(92) 



Generalizing these considerations to a number of k size 
classes, one finds 

= 1 - (1 - (93) 

According to [23 this approach requires a minimal size 
ratio of 7 to 10 between consecutive size classes in order to 
avoid interactions between the packings of the individual 
size classes. 

Discussion of the Furnas model 

The advantages of the FURNAS model are its simple 
derivation, the applicability to an unlimited number of size 
classes and the simplicity of the result. Remarkably, the 
maximum packing fraction does not depend on the volume 
fractions of the individual size classes but only on their 
number. This is the crucial disadvantage of the FURNAS 
model. As an illustration, imagine a given polydisperse 
suspension with a total volume fraction of (^i < (Pxi- In 
the first case, we add a large amount of small spheres, in 
the first case a small amount. The Furnas model would 
yield the same value of the maximum packing fraction 
for both cases, which runs contrary to intuition. In this 
model the maximum packing fraction is calculated without 
effectively considering the particle size distribution in the 
suspension. 



2. An improved model for the maximum packing fraction 

We will derive the new model in a way that at first 
the vivid cases of bi- and tridisperse suspensions are 
treated explicitly in order to prepare the subsequent gen- 
eral derivation. The general formulae will contain the bl- 
and tridisperse suspensions as special cases. 

Description of the volume fractions 

A consequent calculation of the maximum packing frac- 
tion requires the consideration of the individual particle 
size classes' volume fractions. This can be achieved by re- 
taining the particle size distribution existing in the fluent 
suspension state during the composition of the maximum 
packing fraction. This assumption can be expressed in 
the form 



k+l 



(94) 



where 9?^ in the state of maximum packing fraction cor- 
responds to in the fluent state, that is the volume 
fraction of the fcth size class. Introducing the abbreviation 



A4>1 



(95) 



A well as 

4^-1 + 



(96) 
(97) 



Derivation for bidisperse systems 

Bidisperse (or bimodal) suspensions contain two differ- 
ent particle sizes. Given the ratio 



^2 



A0| ^ A0| 



(98) 



we ask how the state of maximum packing fraction can be 
achieved retaining the volume fraction ratio ( [98| . There 
are tow possibilities differing with respect to the value 
of A.(j)\ ■ Both situations are visualized in Figure [s] 












(b) 



Figure 5. Schematic representation of the two possible sit- 
uations for the bidisperse maximum packing fraction: (a) 
Situation 1 - The fraction of large particles lies below ipc while 



the small particles fill the interstices with a fraction of tpc- (b) 
Situation 2 - The large particles occupy a fraction of ipc while 
the small particles are present in the interstices according to 
the volume fraction ratio. 

a. Situation 1 

The volume fraction ip\ of the large spheres lies below 
the monodisperse packing fraction ipc and the remaining 
interstice 1 — Lp'2 is fllled to an amount of ip^. with small 
spheres, compare Figure 5(a)[ So the total particle volume 
fraction is 



VT2 



(99) 



(compare equation (92), where restrictively (f\ — Vc) 
Using the equation ( 98 ) we find from the equation ( 99 1 



^2 



4>l 



(100) 



and furthermore 



fT2 



\A^^ 



(101) 
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In the last step of the equation ( 101 1 we have used the 
identity ( |97|. The hmitation for the validity of the 
equation ( |101[ ) follows from the realizability condition 
</?2 < '^c (arbitrarily, we assign the equal sign to the sit- 
uation 2 in order to avoid an ambiguous definition of 
the case Lp\ = ipc) which in combination with the equa- 
tion yields 



< 



1 



(102) 



as the condition for the validity of situation 1. 
h. Situation 2 

In this situation the volume fraction 1^92 of the large 
spheres equals the monodisperse maximum packing frac- 
tion ipc (it is not possible to exceed the value of ipc in 
a monodisperse loading) and so we have — ipc, see 
Figure |5(b)[ Using the equation ( 98 1 we immediately find 



and 



A(f): 



1 \ V2 

= 



(103) 
(104) 



It follows from the fundamental condition < 1 and 



the equation ( 104 1 that 



1 - V'c 



(105) 



whereas the realizability condition Lp\ < (1 — (^^(/Sc under 
consideration of ip\ = ^Pc and the equation ( 103 I yields 



1 - </5c 



(106) 



The condition (106 1 includes the inequality (1051. So the 



equation ( 106 1 serves as a condition for the validity of 



situation 2. 

c. Summary for bidisperse systems 



classes, there are always two situations. This relies on the 
geometric fact that no particle class except the one with 
the largest particle diameter can reach the monodisperse 
packing fraction. If hypothetically a particle class with 
a smaller diameter reached the monodisperse maximum 
packing fraction, the larger particles would not fit into 
the interstices between the smaller particles and could 
thus not contribute to this state of maximum packing 
fraction. 

a. Situation 1 

The volume fraction if '2 of the large spheres is smaller 
than the monodisperse maximum packing fraction Lpc 
and so the remaining interstice 1 — (f\ \s, filled by the 
medium and smallest particles to an amount of (^^2 (equa- 



tions (101) or (104l) which is the bidisperse maximum 
packing fraction calculated previously. The volume frac- 
tion of the medium and small spheres is thus given by 



^\+^2 = {'^-Vl)V>T2 



(107) 



Rearranging the equation ( 107 1 and using the equa- 
tion (94 1 we find 



(108) 



It is now possible to insert the equation (108 1 into the 
relation ( 107 1 which yields 



^T2 



(109) 



total volume fraction of spheres 



From the equations (|108| and (|109|) we can deduce the 

(110) 



The inequalities (1021 and (106) cover all possible Using equation ( 109 ), we can write the realizability con- 
values of A02. On the basis of a given volume frac- dition if \ < ipc in the form 
tion A02 oi^e has to decide whether situations 1 or 2 



are present. Subsequently the respectively valid equa- 



tions (101 1 or (104) may be used to calculate the maxi- 



ipc(j)2 



mum packing fraction (f^rp^- 

The relations (101) and (104) are — if notation is 
adapted — identical to the equations (3) and (4) presented 
in [B]. 

Derivation for tridisperse systems 

The following considerations refer to tridisperse sys- 
tems, that is systems with three particle size classes. The 
pattern underlying the derivation will later be turned 
into a general formulation for an arbitrary number of size 
classes. 

With tridisperse systems we also distinguish between 
two situations because, regardless of the number of size 



"^.2(1 - V'c) 



(111) 



The inequality (111 I serves as a condition for the validity 
of situation 1 . 
b. Situation 2 

Analogously to the bidisperse case, situation 2 is char- 
acterized by a volume fraction Lp\ = ipc of the largest 
spheres. Under consideration of the equation ( [94| (see 
also equation (IO81) we have 



J® 



fT3 = fi + 'P2 + fa = 



\Aq^l Acj^l 



(112) 



Act^l 
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for the total volume fraction of spheres in situation 2 of a 
tridisperse suspension. The bound for A(/)| follows from 
the condition ip\ + < (1 — Lpc)'f\^2 '^^^ form 



corresponding to inequality ( |111[ ) as expected, 
c. Summary for tridisperse systems 



(113) 



By means of the inequalities ( 111 I and ( 113 ) one decides 



whether the situations 1 or 2 are present. Afterwards, Lplp^ 
is calculated according to the equations (110) or (112)). 



Generalization for polydisperse systems 

Now we will generalize the model for polydisperse sys- 
tems which corresponds to deriving expressions for the 
total volume fraction ^Pt};j^i in the situations 1 and 2 and 
for the limiting value of Ac/))^^-^. 

a. Situation 1 

In the polydisperse case the interstices between the 
largest particles are filled by the smaller particle size 
classes with the packing fraction ff^^ which is determined 
by the volume fractions A<^^ to A(/)^. So the volume 
fraction of the largest particles is V'fe+i and the one of all 
the smaller particles is (1 — V'fe+i)¥'Tfc- sum. of both 
fractions yields the maximum packing fraction, that is 



+ - Vl+l)V^Tk 



(114) 



So we may express the volume fractions of all the small 
particles by 



Inserting the equation ( 118 I into ( 119 I, we find the result 



for the maximum packing fraction in situation 1: 



(120) 



The limiting value for Ac/)^^-^ follows analogously to the 
bi- and tridisperse cases from the requirement < (pc 

and is thus given by 



A0: 



fc+i 



< 



(121) 



b. Situation 2 

In situation 2 the largest particles occupy a volume 
fraction equal to the monodisperse maximum packing 
fraction, so ^p]^_^_l = ipc- Therefore, by means of the 
equation (119) we may write 



•^Tk+l = V>c 



(122) 



for the maximum packing fraction in situation 2. The 
limiting value of Acj)]^^-^ in situation 2 is determined by 
the condition 



(123) 



^ V3„i = (1 - ^pl+iWrk 



m— 1 



Equation ( 94 1 allows for representing all the occuring vol- 



ume fractions by the fraction f^j^i of the largest particles, 
so we can restate the sum in the equation (115) with the 
help of the equation ( 95 1 : 



E 



"Mfc — 



m— 1 



(116) 



It immediately follows from the equations ( 115 ) and ( 116 1 
that 



(1 - '^^k+iWrk = fk+i A ,f 

^'Pk+l 



that can be rearranged into 



VTk^<l>i+l + 't^k 



(117) 



(118) 



Additionally, under consideration of the equation (116) 



the following relation holds for both the situations 1 and 2: 

VTk+l = ^k+l + E ^™ " ( ^ + ) 

"1=1 \ ^k+l J (IIQ'^ 



^k+1 



k+l 



(115) and combined with (115) and (122) takes the form 



A</>: 



fc+i 



> 



(124) 



c. Properties of the bounds 

The bounds (121) and (124) can alternatively be de- 



rived by equaling the prescriptions (120) and (122). This 



underlines the consistency between the situations 1 and 2 
because there is a continuous transition. This is equiva- 
lent to choosing the smallest value out of the two calcu- 
lated maximum packing fractions because for all values 
of Af/))^,^-^ the situation that is present is always the one 
with the smaller maximum packing fraction. Thus we are 
lead to the prescription 



ifiTk+i = min [ffk+i^vfk+i 



(125) 
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F. Scheme for the viscosity estimation 

Figure [6] shows the procedure for the calculation of 
maximum packing fraction and relative viscosity after 
the ith construction step. The complete suspension shall 
contain n particle size classes. 

( start i =1) 



A(t)l from 
equation (74) for 
k = 1, . . . ,i 




Yes 




'Pti 





ip\^2 from equa- 
tion (120) or (122) 



T 



ipxi from equa- 
tion (120) or (122) 




Figure 6. Scheme for the calculation of maximum packing 
fraction and relative viscosity after the ith construction step 
for n particle size classes with references to the respective 
equations 



G. Influence of polydispersity on relative viscosity 

In this section we examine the change in relative vis- 
cosity due to the number of particle size classes, provided 
that the size distribution is given. As already noted in 
section p^VC[ the model developed so far and presented 



in the scheme on page [16] is only valid for diameter ra- 
tios between consecutive size classes of about 7 to 10. A 
special choice of the diameter ratio — in the following we 
will choose Ui = 10 — results in a logarithmic sampling 
of the continuous size distribution. Therefore, the sam- 
pling points, that is the first diameter c?i, have to be 
adequately chosen in order to reproduce the character- 
istics of the distribution. If the volume fractions A(j)k 
are not equally distributed among the different diameter 
values, the choice of the sampling points thus depends on 
the number of size classes. This will be the case in the 
following examination concerning the Rosin-Rammler 
distribution. Subsequently, we will consider the case of a 
uniform distribution. 



1. Rosin-Rammler distribution 

The Rosin-Rammler distribution is frequently used 
to describe the drop size distribution of sprays. Its cumu- 
lative distribution function -Rcdf is given by 



RcBF{d) = 1 - e 



(126) 



where d is the particle diameter as well as X and q are 
model parameters. By differentiation of ( 126[ ) the proba- 
bility density function i?pDF (we do not rigorously distin- 
guish between probability and relative frequency) can be 
deduced: 



RpDF{d) 



(127) 



The function RpDp{d) describes the relative frequency of 
occurrence of a particle diameter between d and d + dd 
and is normalized to unity. So we find the volume fraction 
occupied by the particles with diameters between d and d+ 
dd through multiplication of the total volume fraction 
by i?PDF(rf) (it is 4> = 4>n ^ Z]m=i^'/'™ for ^ number 
of n size classes) . Figure [t] shows the probability density 
function -RpDF (d) for a special choice of the parameters X 
and q in equation ( |127[ ). With a diameter ratio of = 10, 
the distribution only allows for useful sampling at no more 
than three points because of its asymptotic declination. 
The samplings for two and three particle size classes are 
shown in Figure |7] We proceed as follows: In order to 
distribute the total volume fraction (/) among the size 
classes we divide the individual values of the PDF at the 
sampling points by their sum and multiply them by (jj. 
So we achieve the volume fractions to be distributed 
according to the continuous size distribution. Of course, 
it is not necessary to choose sampling points for the 
monodisperse suspension. 

We now evaluate the assignment between the diame- 
ter dk and the relative volume fraction A(f)k given by the 
Rosin-Rammler distribution (127 1 for different values 
of (p. Thereby we apply the e quation (90 I for the viscosity 
as well as the equations ( 120 ) and ( 122 1 for the maximum 



packing fraction. We therefore assume the relation (90 1 
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Q 




100 150 200 250 
particle diameter d in iim 



350 



only causes relatively small changes. The curves diverge 
at the respective values of the final maximum packing 
fraction (pJj^^ which are also depicted in Figure |8] 



2. Uniform distribution 

The uniform size distribution is characterized by the 
equality of all the volume fractions Acjjk, so one has Acjjk — 
(f)/n for each size class. Analogously to the case of the 
Rosin-Rammler distribution, we calculate the relative 
viscosity as a function of the total volume fraction (f) using 
the equation ( |90| and the maximum packing fraction 
applying the relations (1201 and ( |122[ ). The results are 
depicted in Figure |9] whereas the behavior qualitatively 
coincides with the results shown in Figure [8] The relative 



Figure 7. Probability density function by Rosin-Rammler 
(equation ( |127[ )) with parameters X = 100 and g = 1.5 and 
sampling for bi- and tridisperse representations (n — 2: 
[di,d2] = [5,50]pm and n = 3: [di,d2,d3] = [1, 10, 100] pm, 
respectively) 



to be valid and so we de facto apply the FARRIS-Modell 
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(see section IV D 3 1 . 

The results are presented in Figure [8] It is clearly 



300 



250 



9- 



9- 



200 



V 150 
■g 100 



50 




I I 
..' / 
.' / 
.// 
■■ / 



0.2 0.4 0.6 0.8 1 

total volume fraction 

Figure 8. Relative viscosity as a function of the volume fraction 
for various numbers of size classes n from equation | |90[ ) with 
parameters (/Pc = 0.64, \r\\\ = 2.5 and [772] = 6.17; Rosin- 
Rammler size distribution (samplings shown in Figure [t]); 
maximum packing fraction ^'t2 ^^nd (^^3 

visible that the relative viscosity decreases as the number 
of particle size classes increases. The strongest influence 
can be noticed during the transition from a monodisperse 
to a bidisperse suspension while adding a third size class 




0.2 0.4 0.6 O.i 

total volume fraction cj) 

Figure 9. Relative viscosity as a function of the volum e fra ction 
for various numbers of size classes n from equation ( 90 1 with 
parameters ipc = 0.64, [771] = 2.5 und [772] = 6.17 for a uniform 
distribution 

viscosity decreases at constant (j) and increasing n. We 
confine ourselves to at most four size classes because the 
trend is already clearly visible with that number. 



3. Range of influence of polydispersity 

Table |V] shows the viscosity decrease as a function of 
the polydispersity by means of different values of (f) related 
to the monodisperse viscosity. The numerical values for 
the Rosin-Rammler distribution rarely differ from the 
values for the uniform distribution. At a total volume frac- 
tion of = 0.1 the deviation amounts to approximately 
one percentage point, for (j) = 0.2 partially to more than 
five percentage points. If, for instance, the derivation 
shall be kept below one percentage point, it is necessary 
to calculate the viscosity using the polydisperse formulae 
at volume fractions higher than 0.1, provided that the 
particle size distribution is broad enough. A polydisperse 
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Distribution Ratio in % (j) = 0.1 0.2 0.3 0.4 0.5 0.6 



'7r(" = 2) 

ROSIN-RAMMLER Vr(n = l 



100.0 99.1 95.4 87.0 71.1 43.6 7.1 



"•-["-f 100.0 98.7 94.0 83.9 66.1 38.2 5.5 



Uniform Distribution nrlfl^ iqq.O 98.6 93.6 83.4 65.7 38.6 5.9 



100.0 99.0 95.1 86.9 71.9 46.2 8.6 
100.0 98.6 93.6 83.4 65.7 38.6 5.9 
100.0 98.5 92.9 81.7 63.0 35.5 5.1 



Table V. Ratio of the relative viscosities of polydisperse {r]r{n > 2))and monodisperse systems (rjrin = 1)) with given size 
distribution, data in percentage points related to the monodisperse viscosity for various volume fractions, calculated from 



equation 1 90 1; parameters: tpc ~ 0.64, [»7i] = 2.5, [772] ~ 6.17 



calculation only makes sense if the size distribution allows 
for at least two size classes with significant volume frac- 
tions having a diameter ratio of 7 to 10. For example, the 
inclusion of a fourth size class with a particle diameter 
of ^4 = 1000 pm in the case of the Rosin-Rammler dis- 
tribution in Figure [7] has no effect on the relative viscosity 
because the relative frequency is approximately equal 
to 9 X 10"^^ and therefore negligible. 



The equation (130 1 allows for the effect of the particle 



viscosity r]p on the relative viscosity of polydisperse sys- 
tems to be described approximately. Since in this work 
we do not focus on particle deformation, we simply state 
the result ( 130 1 without further validation. 



CONCLUSIONS 



H. Consideration of particle deformation 

In [21] an approach to the consideration of deformable 
particles is developed. Thereby, the particle deformability 
is represented by a particle viscosity rjp. Using the so- 
called modified ESHELBY model (see also |H]) as well as 
the elastic- viscous analogy the derivation of the equation 



Vr 



2?7p + 3?7o + 30 {r]p - rjo) 
2r]p + 3t]o - 2(j) {rip - 770) 



(128) 



is outlined. This equation is confined to long-range parti- 
cle interactions and is thus invalid for dense suspensions. 
However, as presented in [2T|, the first-order Taylor 
series expansion of the relation ( |128[ ) can be used for the 
differential Bruggeman model (see section IV A I, analo- 
gously to the use of the Einstein relation in equation ( 34 1 . 
So for the differential approach we have 



r] + dr] = 1] 



l + 2.5(/./ 



7]p + 1.577 



(129) 



To be consistent with the equation ( |128 l, one has to 
set [7^ 1] = 2.5. For 77p — >■ 00 (rigid particles) the equa- 
tion (1291 reduces to the relation ([Ml. In the course 



of the discrete construction process higher order terms 
out of the series expansion of the equation (128) can be 



considered. However, the limited validity of this equation 
for larger values of 0p must be regarded. If we include 
terms up to the second order we can express the relative 
viscosity in the discrete construction process by 



1 -t- 2.5(/ip,i+i 



2-5 0p_,+i 



■rip - iq 

TJp + 1.57/j 



rjp + 1.577,; 



(130) 



In the present work we provided a model for the rel- 
ative viscosity of polydisperse suspensions of spherical 
non-colloidal particles. Using monodisperse viscosity cor- 
relations, we described polydisperse suspensions by means 
of a construction process consisting of successive additions 
of particle size classes. 

As a starting point, we proposed a generalized form of the 
well-known Krieger-Dougherty equation that allows 
for the choice of the second order intrinsic viscosity [772]- 
This modified equation can be used to approximate the 
various monodisperse viscosity relations existing in the 
literature and can therefore be regarded as a generic rela- 
tion. 

Later, we described the construction process in detail 
applying a dimensionless way of description based on 
volume fractions. This rigorous description served as a 
basis for the calculation of the relative viscosity during 
the construction process. Starting from the Bruggeman 
model, we finally arrived at the Farris model, connecting 
two approaches commonly regarded as uncorrelated. 
As an entirely new component, we introduced the polydis- 
perse maximum packing fraction into the Farris model. 
Here, we followed a physically consistent approach in 
contrast to the approaches presented in the literature. 
Consistently with the relative viscosity calculation, we 
derived a formalism to determine the polydisperse maxi- 
mum packing fraction by means of a common construction 
process. 

The entire formalism for calculating the relative viscosity — 
whether including the maximum packing fraction or not — 
is depicted on page [16] We evaluated the model in the 
case of two different particle size distributions in order to 
observe the influence of polydispersity represented by the 
number of particle size classes. 

Additionally, we revealed a possible approach for integral- 
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ing particle deformability, represented by a particle phase 
viscosity, into the viscosity model using a result from the 
literature. 

So far, our model is only valid for large diameter ratios of 
consecutive size classes during the construction process. 
An attempt to generalize the model to the case of small 
size ratio as well as the shear rate dependence of the 
relative viscosity will be presented in a future work. 



is given by 



(A3) 



while the volume fraction (/)^°"^' from the equation ( 75 1 is 



1 - Em=i+1 



(A4) 



Appendix A: Proof of equivalence of the 
construction processes at variable and constant 
volume 



The Definition (46 1 implies for both cases — variable 



and constant total volume — the corollary 



Insertion of the relation (|60| into the equation ( A4 1 yields 
after an intermediate step 



(A5) 



E 



(Al) 



and so the validity of the proposition ( A2 1 has been proved 



It follows directly the coincidence of the differences S(j)i+i 
in the equation (46 1. In a second and last step we proof 
that 



Obviously both approaches must result in the same total 
volume fraction at the end of the construction process. 
So we ask if all partial sums 4>i {i < n) are identical in 
both cases. As a consequence, the values of S(pi+i for 
i = 1 . . .n would coincide, too. The validity of these 
identities is important for the derivations in section [IV D| 
because it ensures that the viscosities calculated in both 
cases are equal. This is a prerequisite for a consistent 
model. In the following we show the equivalence of the 
two approaches. 



To distinguish between the volume fractions 



oc- 



curring in each case we introduce the notations (1)*^^ for 
the case of variable total volume and ^i*^™*'* for the case 
of constant total volume. In a first step we show that the 
volume fraction (/), coincide for i = 1 . . .n, so that 



(A2) 



/ var 



I* var 



1 + ^tTi 



■'p,i+i 



(A6) 



of the equation ( A6 ) in the form 



Using the definition ( 43 1 we may write the left-hand side 



x*var 



(A7) 



Rearranging of the right-hand side of the representa- 
tion ( A6 ) and using the equations ( 70 1 and ( 60 1 we find 



i^const 



(A8) 



var 



Vf + E 



4 + 1 



1 + ct>nr 



According to the equation ( 42 1 , the volume fraction 



So the proof of the proposition ( A6 1 is complete 
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